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ABSTRACT 

We report our discovery of orbitally modulated y-ray emission from the black widow system PSR 
J13II—3430. We analyze the Fermi Large Area Telescope data during the offpulse phase interval 
of the pulsar, and find the orbital modulation signal at a ^3 (t confidence level. Further spectral 
analysis shows no significant differences for the spectra obtained during the bright and faint orbital 
phase ranges. A simple sinusoid-like function can describe the modulation. Given these properties, 
we suggest that the intrabinary y-ray emission arises from the region close to the companion and the 
modulation is caused by the occultation of the emitting region by the companion, similar to that is 
seen in the transitional millisecond pulsar binary (MSP) PSR J1023-I-0038. Considering the X-ray 
detection of intrabinary shock emission from eclipsing MSP binaries recently reported, this discovery 
further suggests the general existence of intrabinary y-ray emission from them. 

Subject headings: binaries: close — pulsars: individual (PSR J1311—3430) — gamma rays: stars 


1. INTRODUCTION 

Millisecond pulsars (MSPs) are widely accepted to 
be old neutron stars that were spun up through 
mass accretion from the companions whe n they were 
at t h e low-mass X-ray binary phase (jAlnar et al.l 
119821 : iRadhakrishnan fc Srinivasanl Il982fl . Not sur- 
prisingly, >60% know n MSPs are in binaries (e.g., 
[Manchester et al.ll2005h . A s ub-class of them, so -called 
‘black widow’ pulsar systems (|Fruchter et ahlfTOSSll , have 
very low-mass, ^0.02 Mq companions. To form isolated 
MSPs, one possible channel is through ablation of the 
companions by the pulsar wind. This possibility likely 
occurs in the black widows since they are eclipsing sys¬ 
tems at radio frequencies, indicating the interaction be¬ 
tween the pulsar wind and the companions. X-ray ob¬ 
servations of them revealed orbital flux variations, also 
supporting the presence of the intrab inary interaction 
(|Huang et al.ll20l3lGentile et al.ll2014ll . In addition, re- 
cent extensive studies of the so-called ‘redback’ systems 
(|Robertsll2013ll have provided clear evidence for the in¬ 
teraction. These redbacks are also eclipsing MSP bi¬ 
naries, but contain relative massive, '^O.I-O.GMq com¬ 
panions. X-ray observations of the prototypical red- 
back PSR J1023-I-0038 detected significant orbit al flux 
variations (lArchibald et al.l[2010HBogdanov et al.ll2011ll . 
and the variations can be expl ained by the existenc e 
of an intrabinary shock region (iBogdanov et al.l 1201111 . 
Similar features w ere also clearly seen i n the redback 
XSS J12270—4859 (|Bogdanov et al.ir2014l and references 
therein). 

Owing to its all-sky monitoring and high sensitiv¬ 
ity capabilities, the Fermi Gamma-ray Space Telescope, 
launched in 2008, has greatly improved our studies of 
pulsars. For MSPs, more than six-fold black widows and 
redba cks have been discovered with the help of Fermi 
(e.g.. IRobertsll201^ . At Fermi's lOOMeV to 300GeV 
energy range, marginal evidence for the intrabinary in¬ 
teraction in the eclipsing systems has also been seen. 
For the first discovered black widow PSR B1957-I-20 


(|Fruchter et al.l[198811 . an orbital modu lation signal wa s 
detected at a ^ 2.3a conhdence level (|Wu et al.l 1201211 . 
In addition, p ossible signals were also reported for XSS 
J12270—4859 (|Xing fc Wang|l2014ll and a can didate red- 
back 2FGL J0523.3-2530 Jxing et aLllMl . Theoreti¬ 
cal studies have long predicted the intrabinary interac¬ 
tion and related high-energy emission from black widows 
(e.g. JArons fc Tavanil 119931 . Studies of the 7-ray emis¬ 
sion from the intrabinary region allow us to explore the 
detailed physical pr ocesses within such a binary (e.g., 
[Roberts et al.l[2M3l . In this paper we report the de¬ 
tection of orbitally modulated 7-ray emission from a re¬ 
cently discovered black widow PSR J1311—3430, which 
thus indicates the intrabinary origin for part of its emis¬ 


sion. 


PSR J1311—3430 was initially listed as an unassoci¬ 
ated source in the Fermi Large Are a Telescope (LAT ) 
source catalog (2FGL J1311.7— 3429: lNolan et al]l2012ll . 
It is the only 7-ray selected MSP with 7-ray pulsed emis- 
sion discovered via a direct blind search in the Fermi data 
(|Pletsch et al.ll2012f). The pul sed radio emission soon was 
detected too (|Rav et al.ll2013ll . but the signal was visible 
only during <10% of the observation time, suggesting 
strong variations in the intrabinary medium. Before the 
discovery, the source was found to have orbital modula¬ 
tion with a short period of :s:94 minutes through optical 
imag ing and spectroscopy (|Romanil 120121 : IRomani et al.l 
1201211 . Gonsidering its properties of weak X-ray emis¬ 
sion, sinusoid-like optical modulation and large modula¬ 
tion amplitude, it was already suggested to be a black 
widow system and the optical modulation is c aused by 
irrad iation of the companion by the pulsar wind (|Romanil 
nnn). Marginal modulated X-ray emission possibly re- 
lated to the intrabinary shock has also been detected 
(|Romanil 120121 : iKataoka et al.l 1201211 . The 7-ray discov¬ 
ery of the 2.5 6 ms spin signal th us confirmed its black 
widow nature (|Pletsch et al.l[2M^ . Analyzing the Fermi 
data, we searched and found the orbital modulation sig¬ 
nal from the source’s offpulse emission. Below we present 
the data analysis and results in Section O The results 
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are discussed in Section [31 
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Fig. 1.— Folded pulse profile and two-dimensional phaseogram 
in 32 phase bins obtained for PSR J1311—3430. The gray scale 
represents the number of photons in each bin, and the dashed lines 
mark the onpulse and offpulse phase intervals. 


2. DATA ANALYSIS AND RESULTS 

2 . 1 . Fermi LAT data 

We selected 0.1-300 GeV LAT events from the Fermi 
Pass 7 Reprocessed (P7REP) database inside a 20° x 20° 
region centered at the position of PSR J1311—3430 
(jPletsch et al.l[2(TT^ during the time period from 2008- 
08-04 15:43:36 to 2015-01-06 21:19:57 (UTC). Only 
events with zenith angle less than 100 deg and during 
good time intervals were kept. The former prevents the 
Earth’s limb contamination, and for the latter, the qual¬ 
ity of the data was not affected by the spacecraft events. 


2 . 2 . Timing Analysis 

We performed timing analysis to the 0.1-300 GeV LAT 
data of the PSR J131 1—3430 region to up date the 7 - 
ray ephemeris given in iPletsch et al.l (l2012f) . An aper¬ 
ture radius of 1?0 was used. We determined the pulse 
time of arrivals (TOAs) by obtaining the pulse profiles of 
40 evenly divided se gments using the known ephemeris 
(jPletsch et al.l l2012f l and cross-correlated them with a 
template profile created with data during the time pe- 
ri od of MJD 54682-5 6119 (the same time range as that 
iniPletgch et^jl.| [2012ll . following the a lgorithm described 
iniTavloil (1199211. We used TEMP02 (iHobbs et al.l[200a 
lEdwards et al.l 200611 to fit the TOAs. Only pulse fre¬ 


quency / and frequency derivative / were fitted, and 
the other timing parameters were fixed to their known 
values. We obtained / = 390.56839326403(7) Hz and 
/ = —3.193(1) X 10“^® s“^, consistent with the values 


given in IPletsch et al.l (| 2012 ll within ^O.Scr and ^ 2 . 2 cr 
uncertainties, respectively. The folded pulse profile and 
two-dimensional phaseogram are shown in Figure [TJ We 
defined phase 0.16-0.66 and 0.66-1.16 as the onpulse and 
offpulse phase intervals, respectively. 



Orbital Phase 


Fig. 2.— 0.2—300 GeV light curve folded at the orbital period 
using the offpulse data. The bottom (Phase I) and peak (Phase 
II) ranges of the modulation are marked. Two simple sinusoid fits 
are displayed as dashed and dotted curves; for the latter, 0.1 phase 
shift is forced (see the text in Discussion). 

2.3. Maximum Likelihood Analysis 

We selected LAT events in 0.1-300 GeV energy range 
for the likelihood analysis, and included all sources 
within 20 deg in the Fermi thi rd source catalog 
(|The Fermi-LAT Collaborationl20'T5l) centered at the po¬ 
sition of PSR J1311—3430 to make the source model. 
The spectral function forms of the sources are provided 
in the catalog. The spectral parameters of the sources 
within 5 deg from PSR J1311—3430 were set free, and all 
other parameters of the sources were fixed at their cat¬ 
alog values. The 7 -ray counterpart of PSR J1311—3430 
was modeled with an ex ponentially cutoff power law, 
characteristic for pulsars (lAbdo et al.ll2013ll . and a sim¬ 
ple power law for comparison. In addition, we used the 
spectrum model glLiem_v05_revl.fits and the spectrum 
file iso_source_v05.txt to consider the Galactic and ex- 
tragalactic diffuse emission, respectively. 

Using the LAT science tools software package v9r33p0, 
we performed standard binned likelihood analysis to the 
LAT data. The 7 -ray emission during the total pulse 
phase interval was detected with a Test Statistic (TS) 
value of 6279, while that during the onpulse and off¬ 
pulse phase intervals were detected with TS values of 
7723 and 499, respectively. The TS value at a specific 
position is calculated from TS= —2 log(Lo/ii)) where Lq 
and Li are the maximum likelihood values for a model 
without and with an additional source respectively, and 
approximately is the squa re of the detection significance 
for the additional source (jAbdo et al.ll20inD . We found 
during the total pulse phase, onpulse phase, and offpulse 
phase intervals, the emission is better modeled by an ex¬ 
ponentially cutoff power law, with the low energy cutoff 
detected with >13cr, >14cr, and >5a significance (esti- 
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Fig. 3.— 0.2—300 GeV TS maps of a 2° x 2° region centered at the position of PSR J1311—3430 during Phase I (left) and Phase II 
[right). The image scale of the maps is 0.04° pixel”^. The color bars indicate the TS value range. All sources in the source model except 
PSR J1311—3430 were considered and removed. The white [left) and dark [right) crosses mark the position of PSR J1311—3430, and the 
green cross marks the position of the nearby catalog source 3FGL J1316.0—3338. 


mated from y—2 log(Lp;/Lexp), where and Lpi are 
the maximum likelihood values for the exponentially cut- 
off power-law mo del and power-law model, respectively; 
lAbdo et all l2013fl . The resulting exponentially cutoff 
power-law fits are summarized in Table [T] 

2 . 4 . Orbital Variability 

We folded the LAT events of the PSR J1311—3430 re¬ 
gion at its orbital period (iPletsch et al.ll201^ to study 
its poss ible orbital rnodula tions. The source position 
given in iPletsch et al.l (|2012t) was used for the barycentric 
corrections to photon arrival times, and photons within 
Rmax (Rmax ranges from 0?1-1?0 with a step of 0?1) from 
the position were collected. Different energy ranges (0.1- 
300, 0.2-300, 0.3-300, 0.5-300,1-300 GeV) were tested in 
folding. No significant modulations were detected using 
the whole data (i.e., during the total pulse phase), as the 
largest H-test val ue was 6 (correspond ing to <2cr detec¬ 
tion significance: Ide Jaeer et al1ll989ll . However, a sig¬ 
nificant orbital signal was best revealed using the offpulse 
data in the >0.2 GeV energy range within 0?4 from PSR 
J1311—3430. The folded light curve, which has an H- 
test value of ~22 (corresponding to ^4cr significance and 
^3(7 post-trial significance, for the latter where 50 trials 
on the energy range and aperture radius are considered), 
is shown in Figure [5] The phase zero is a t the ascend¬ 
ing node of the pulsar (|Pletsch et al.ll20lll . The folded 
light curve has a brightness peak around the superior 
conjunction (when the companion is behind the pulsar), 
the sa me as the modest X-ray one reported in IRomanil 
(|2012f) . The similarity helps strengthen the y-ray mod¬ 
ulation detection. Using the LAT tool gtexposure, we 
checked t he summed exposures over the 10 orbital phase 
bins fe.g.. [Johnson et ^l2015[l . and they had only <1% 
differences, too small to cause any artificial orbital mod¬ 
ulations. 

We performed likelihood analysis to the >0.2 GeV off¬ 
pulse LAT events during the orbital phase ranges of 0.2- 
0.5 (named Phase I) and 0.7-1.0 (named Phase II), which 
were approximately defined for the bottom and peak of 
the orbital modulation, respectively. We found that the 


emission during the both phase ranges is better modeled 
by an exponentially cutoff power law, with low energy 
cutoff detected with >3cr and >2 (t significance. The ex¬ 
ponentially cutoff power-law fits are summarized in Ta¬ 
bled] The TS values during Phase I and II are ~120 and 
^320 (see Figure[3|), respectively, which indicate that the 
source during the latter is more significantly detected 
than during the former, confirming the detection of or¬ 
bital modulation from photon folding. 

Possible contamination from a nearby catalog source 
3FGL J1316.0—3338, which is identified as the coun¬ 
terpart to the f lat spectrum radio q uasar (FSRQ) 

PKS 1313-333 (|Ackeman et al.l was investi¬ 

gated. This source, being only ^1?2 away from PSR 
J1311—3430 and relatively bright (TS~625 in the cat¬ 
alog), exhibited flaring events in the pasiQ. Using the 
offpulse data and performing likelihood analysis, we ex¬ 
tracted its 30-day interval light curve, and found that for 
five time bins (MJD 54863-54923, MJD 55433-55493, 

MJD 56363-56393), it had fluxes >2 ct above the value 
obtained from the total offpulse data. We repeated the 
analysis by excluding the data of the time bins. We found 
that the folded light curve is nearly the same, still having 
an H-test value of ~22. 

2 . 5 . Spectral Analysis 

We further investigated the orbital-dependent spectral 
variability during the offpulse phase interval. Spectra 
of PSR J1311—3430 during the whole offpulse phase in¬ 
terval, Phase I, and Phase H were obtained, and the 
spectrum during the onpulse phase interval was also ob¬ 
tained for comparison. We extracted the spectra by per¬ 
forming maximum likelihood analysis to the LAT data in 
10 evenly divided energy bins in logarithm from 0.1-300 
GeV, with the emission of the source being modeled with 
a power law in each energy bin. We only kept spectral 
points with TS>4, and derived the 95% upper limits in 
other energy bins. The spectra extracted by this method 

^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/4yr_catalog/ap_lcs.php 
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Fig. 4.— Left panel: 7 -ray spectra of PSR J1311—3430 obtained with the onpulse (green triangles) and offpulse data (dark circles). Right 
panel: 7 -ray spectra of PSR J1311—3430 obtained during offpulse Phase I (blue diamonds) and Phase II (red squares). The exponentially 
cutoff power-law fits for the onpulse and offpulse data are shown as green solid and dark dot-dashed curves, respectively. The same model 
fits for the offpulse Phase I and Phase II data are shown as blue dotted and red dashed curves, respectively. 


are less model-dependent and provide a detailed descrip¬ 
tion of the 7 -ray emission for the source. 

The obtained spectra are shown in Figure ID The on¬ 
pulse emission appears to have a 3-times higher cutoff 
energy than the offpulse one (see also Table [IJ . In 
addition, comparing the two offpulse spectra, the source 
was brighter across the >0.2 GeV energy range during 
Phase II than during Phase I. 

We also repeated the analysis by excluding the data 
when the nearby source 3FGL J1316.0—3338 had possi¬ 
ble flares (see §[I 1 - The obtained spectral parameters 
of PSR J1311—3430 during the two orbital-phase ranges 
are consistent with the values obtained above (within Icr 
uncertainties). We concluded that the flares do not have 
any significant effect on our spectral analysis. However, 
we note that in the lowest energy bin, the two upper lim¬ 
its of the Phase I and II spectra (see the right panel of 
Figure S]) are lower than the exponentially cutoff power- 
law fits. This problem might be due to possible con¬ 
tamination from 3FGL J1316.0—3338 in the low-energy 
range. 

3. DISCUSSION 

From our analysis of the Fermi data of PSR 
J1311—3430, we have detected its y-ray orbital modu¬ 
lation during the offpulse phase interval of the pulsar, 
where the magnetospheric emission from the pulsar was 
likely effectively removed. In bo th optical and X-ray ob¬ 
servations, flares were detected (lRomanill2012D . indicat¬ 
ing the strong interaction between the pulsar wind and 
the companion. Likely y-rays are also produced due to 
the intrabinary interaction. However, different from that 
is seen in PSR B1957-I-20, which has an extra component 
above 2.7 GeV at its inferior conjunction (when the com¬ 
panion is in front of the pulsar), the light curve peak 
is near the superior conjunction for PSR J1311—3430. 
The difference suggests that the intrabinary y-ray emis¬ 
sion model, which explains the extra component as the 
result from viewing an in verse Gompton p rocess as a 
he ad-on coll i sion ( see also iBednareU (120141) 1. proposed 
in IWu et'aD (| 2012 D does not apply here. 

For PSR J1311—3430, no significant spectral changes 
were found from the offpulse orbital-phase-resolved spec¬ 


tra (Figure SI). The source appeared brighter across the 
>0.2 GeV energy range during the peak range (Phase H) 
than during the bottom range (Phase I). Although the 
uncertainties from our fits with the exponentially cut¬ 
off power law are relatively large, the two spectra are 
generally similar to each other. The similarity suggests 
an geometric origin for the orbital modulation, such as 
that used to e xplain the X-ray orbita l modulation of PSR 
J1023-I-0038 (iBogdanov et ai.ir201lD . The binary likely 
has an inclination angle of i ^60°, and the companion 
has a very small Roche lobe radius Rl = 0.068 Rq (a 
canonical neutron star mass Mn = 1.35 Mq is assumed). 
Therefore like in PSR J1023-I-0038, the intrabinary emis¬ 
sion region must be very close to the companion and the 
companion can thus block part of the region, causing the 
observed orbital modulation. Here we simply assume 
that the region is at the inner surface of the companion, 
and use a function m/j[l -|- sin(27r^ — tt) sinz ]/2 -|- rric to 
describe the orbital modulation, where (j) is the orbital 
phase, i = 60° is fixed, and mh and TOc are the modula¬ 
tion amplitude (in units of counts) and constant counts, 
respectively. Fitting the folded light curve, we found 
this function can describe the modulation (see Figure [2]), 
where the minimum = 9.7 (for 8 degrees of freedom) 
and rrih = 27 ± 7, TOc = 27 ± 3. However, examining 
the light curve, the minimum and maximum may have a 
0.1 phase shift (i.e., they occur at phase 0.35 and 0.85, 
respectively). If we force such a shift, the results are 
= 7.9 (for 8 degrees of freedom) and ruh = 28 ± 7, 
rric = 26 ± 3, indeed slightly better. This shift may sug¬ 
gest that the emission is not isotropic. The constant part 
rric may represent the intrabinary emission unblocked by 
the companion over the whole orbital phase, while emis¬ 
sion from the pulsar during its offpulse phase interval 
could also contribute a small fraction. 

It is not clear why PSR J1311—3430 has an orbital 
modulation different from that of PSR B1957-I-20. We 
note that if the size of the interaction region is propor¬ 
tional to that of a companion, similar fractions of an 
isotropic pulsar wind (0.0021 vs. 0.0024) would be in¬ 
tercepted for PSR J1311-3430 and PSR B1957-f20 re¬ 
spectively (estimated from (i? 2 / 73 b)^/ 4 ; R 2 is the radius 
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TABLE 1 

Exponentially cutoff power-law fits for PSR J1311—3430. 


Data set 

>0.1 CeV Flux 

(10“® photon cm“^ 

F 

(GeV) 

TS 

total data 

8.0 ± 0.4 

1.80 ± 0.04 

4.0 ± 0.5 

WI7W 

Onpulse data 

13.9 ± 0.6 

1.71 ± 0.04 

3.9 ± 0.4 

7723“ 

Otlpulse data 

4.6 ± 0.6 

1.9 ± 0.2 

1.3 ± 0.4 

499 

Otlpulse Phase 1 

3 ± 1 

1.2 ± 0.7 

0.6 ± 0.3 

mw 

Otlpulse Phase 11 

7T2 

1.8 ± 0.5 

1.2 ± 0.9 

320 


Note. — Column 3 and 4 list the photon index and cutoff energy of the exponentially cutoff power-law model. 


of the companion and Db is the separation distance of 
the binary). The notable differences are that the spin- 
down luminosity of PSR J1311—3430 is approximately 
1/3 of PSR B1957-I-20 and that the companion in PSR 
J131 1—3430 nearly fills its Roche lobe (|Romani et al.l 
1201211 . We suspect that because of the latter, an out¬ 
flow from the companion may still exist, the same a s 
that in PSR J1023-I-0038 (e.g., iBogdanov et al.ir2nilh . 
For PS R B1957-I-20, its cornp anion fills its Roche lobe 
^ 85% l|Revnolds et alJ 1200711 . and the mass loss from 
the companion is presu mably driven by the pulsar wind 
(|Arons fc Tavanill993l and references therein). The much 
more messy environment in PS R J1311—3430, as sug¬ 
gested by the radio observations (jRav et al.ll2013[l , could 
be evidence for this possibility. Detailed modeling for 
physical processes in PSR J1311—3430 would help verify 
it. 

Our discovery of the orbital 7 -ray modulation in PSR 
J1311—3430 and the analysis results have provided clear 
evidence for 7 -ray production due to intrabinary inter¬ 
action between a pulsar and its companion in a black 
widow system, and thus have confirmed the general phys¬ 


ical picture that has long been proposed theoretically. 
Since X-ray observations have revealed the general ex¬ 
istence of intrabinary shock emission in eclipsing MSP 
binaries, similar work can be carried out to search and 
study related 7 -ray emission from these recently identi- 
hed systems. 
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